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Abstract: Renormalization of Hamiltonian field theory is usually a rather painful 
algebraic or numerical exercise. By combining a method based on the many-body 
coupled cluster method, analyzed in detail by Suzuki and Okamoto, with a Wilsonian 
approach to renormalization, we show that a powerful and elegant method exist to solve 
such problems. The method is in principle non-perturbative, and is not necessarily 
unitary. 
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1. Introduction 

The power of Hamiltonian methods is well known from the study of non-relativistic 
many-body systems and from strongly-interacting few particle systems, even though a 
Lagrangian approach is usually chosen for relativistic theories. Hamiltonian methods 
for strongly-interacting systems are intrinsically non-perturbative and usually contain 
a Tamm-Dancoff type approximation, in the sense that, in practice, one has to limit 
the bound state as an expansion over states containing a small number of particles. 
This truncation of the Fock space gives rise to a new class of non-perturbative diver- 
gences, since the truncation docs not allow us to take into account all diagrams for 
any given order in perturbation theory. Therefore rcnormalization issues have to be 
considered carefully. Two very different remedies for this issue are the use of light- 
front Tamm-Dancoff field theory (LFFT) [1] and the application of the coupled cluster 
method (CCM) [2, 3]. In the LFFT the quantization plane is chosen to coincide with 
the hght front, therefore the divergences that plagued the original theories seem to dis- 
appear [4] . Furthermore, not having to include interactions in boost operators allows a 
renormalizable truncation scheme [5]. One of the most important difficulties in LFFT 
is the complicated structure of the rcnormalization process [6]. In the standard form 
of CCM, on the other hand, the amplitudes obey a system of coupled non-linear equa- 
tions which contain some ill-defined terms because of ultraviolet divergences. It has 
been shown [7] that the ill-defined amplitudes, which are also called critical topologies, 
can be systematically removed, by exploiting the linked-cluster property of the ground 
state. This can be done by introducing a mapping which transfers them into a finite 
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representations without making any approximation such as a couphng expansion. Thus 
far this resummation method has been restricted to superrenormahzable theories due 
to its complexity. 

A natural question in the renormalization within the Hamiltonian formalism arises, 
since one could also perform renormalization in a Lagrangian framework and finally 
construct the corresponding Hamiltonian by means of standard Legendre transforma- 
tion. [It should be noticed that this is not generally applicable if the Lagrangian contains 
higher-order time- derivatives.] The Hamiltonian formalism, despite a certain lack of el- 
egance, has the advantage that it is very economical, and one can use all the know-how 
of quantum many-body theory. 

In the last decade extensive attempts have been made to give a prescription for 
renormalization within the Hamiltonian formalism [8, 9, 10]. Commonly unitary trans- 
formations are used to decouple the high- and low-energy modes aiming at the partial 
diagonalization of the Hamiltonian. One of the elegant approach in this context has 
been the so-called similarity renormalization group (SRG) proposed by Glazek and 
Wilson [8] (and by Wegner [9] independently). The SRG resembles the original Wilso- 
nian renormalization group formalism [11], since a transformation that explicitly runs 
the cutoff is developed. In an early paper, Wilson [11] exploited a similar transfor- 
mation which was originally introduced by Bloch [12]. There the Hamiltonian matrix 
is reduced by using a transformation which lowers the cutoff initially imposed on the 
individual states. Later, Wilson abandoned this formalism in favour of a Lagrangian 
one. The most important reason was that the Bloch transformation is ill-defined and 
produces fake divergences. These divergences emerge from denominators which contain 
a small energy difference between states retained and states removed by the transfor- 
mation. The SRG [8, 9] is designed to be free of such small energy denominators and to 
eliminate interactions which change the energies of the unperturbed states by a large 
amount. However, there are several issues in this approach: it is hard to incorporate 
loop expansions within the method, the SRG can not systematically remove interac- 
tions which change the number of particles (i.e, when the Hamiltonian is not diagonal 
in particle number space), and most importantly, the computations are complex and 
there is no an efficient non-perturbative calculable scheme. 

In this paper we introduce a method for obtaining the low energy effective operators 
in the framework of a COM approach. The transformation constructed avoids the small 
denominators that plague old-fashioned perturbation theory. Neither perturbation the- 
ory nor unitarity of the transformation are essential for this method. The method is 
non-perturbative, since there is no expansion in the coupling constant; nonetheless, 
the CCM can be conceived as a topological expansion in the number of correlated 
excitations. We show that introducing a double similarity transformation using linked- 
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cluster amplitudes will simplify the partial diagonalization underlying renormalization 
in Hamiltonian approaches. However, a price must be paid: due to the truncation the 
similarity transformations are not unitary, and accordingly the hermiticity of the re- 
sultant effective Hamiltonian is not manifest. This is related to the fact that we have a 
biorthogonal representation of the given many-body problem. There is a long tradition 
of such approaches. The first we are aware of are Dyson-type bosonization schemes 
[13]. [Here one chooses to map the generators of a Lie algebra, such that the raising 
generators have a particularly simple representation.] The space of states is mapped 
onto a larger space where the physically realizable states are obtained by constrained 
dynamics. This is closely related to CCM formalism, where the extended phase space 
is a complex manifold, the physical subspace constraint function has been shown to be 
of second class and the physical shell itself was found to be a Kahler manifold [14]. The 
second is the Suzuki- Lee method in the nuclear many-body (NMT) problem [15, 16], 
which reduces the full many-body problem to a problem in a small configuration space 
and introduces a related effective interaction. The effective interaction is naturally un- 
derstood as the result of certain transformations which decouple a model space from 
its complement. As is well know in the theory of effective interactions, unitarity of 
the transformation used for decoupling or diagonalization is not necessary. Actually, 
the advantage of a non-unitarity approach is that it can give a very simple description 
for both diagonalization and the ground state. This has been discussed by many au- 
thors [17] and, although it might lead to a non-hermitian effective Hamiltonian, it has 
been shown that hermiticity can be recovered [14, 18]. Nevertheless non-hermiticity is 
negligible if the model space and its complement are not strongly correlated [19, 20]. 
Therefore defining a good model space can in principle control the accuracy of CCM. 

To solve the relativistic bound state problem one needs to systematically and si- 
multaneously decouple 1) the high-energy from low-energy modes and 2) the many- 
from the few-particle states. We emphasize in this paper that CCM can in principle 
be an adequate method to attack both these requirements. Our hope is to fully utilize 
Wilsonian Exact renormalization group [22] within the CCM formalism. Here the high 
energy modes will be integrated out leading to a modified low-energy Hamiltonian in 
an effective many-body space. Notice that our formulation does not depend on the 
form of dynamics and can be used for any quantization scheme, e.g., equal time or 
hght-cone. 

The organization of this paper is as follows. In section II, we discuss our approach 
and it's foundation. Finally we conclude and present an outlook in section HI. 
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2. Formalism 



The discussion in this section is partiaUy based on the work of Suzuki and Okamoto 
[16]. Let us consider a system described by a Hamihonian H{A) which has, at the very 
beginning, a large cut-off A. We assume that the renormahzed Hamihonian H'^^{A) up 
to scale A can be written as the sum of the canonical Hamiltonian and a "counterterm" 

H'''{A) = H{A) + Hc{A). (2.1) 
Our aim is to construct the renormalized Hamiltonian by obtaining this counterterm. 
Now imagine that we restrict the Hamiltonian to a lower energy scale (n), where we 
want to find an effective Hamiltonian H^^{^) which has the same energy spectrum 
as the original Hamiltonian in the smaller space. Formally, we wish to transform the 
Hamiltonian to a new basis, where the medium-energy modes /i < k < A, decouple from 
the low-energy ones, while the low-energy spectrum remains unchanged. We define 
two subspaces, the intermediate-energy space Q containing modes with fi < k < A 
and a low-energy space P with fi < k. Our renormalization approach is based on 
decoupling of the complement space Q from the model space P. Thereby the decoupling 
transformation generates a new effective interaction SH{iJ,, A) containing the physics 
between the scale A and /i. One can then determine the counterterm by requiring 
coupling coherence [10, 23], namely that the transformed Hamiltonian has the same 
form as the original one but with A every where replaced by /i. [This is in contrast 
to the popular Effective Field Theory approach, where one includes all permissible 
couplings of a given order and fixes them by requiring observable computed be both 
cutoff-independent and Lorentz covariant.] The operators P and Q which project a 
state onto the model space and its complement, satisfy P^ = P, = Q, PQ = and 
P + Q = 1. We introduce an isometry operator G which maps states in the P- onto 
the Q- space, 

\q)^G\p) (|g)eg,|p)eP). (2.2) 

The operator G is the basic ingredient in a family of "integrating-out operators" , which 
passes information about the correlations of the high energy modes to the low-energy 
space. The operator G obeys G = QGP, = 0, PC = and = forn ^ 2. The 
rather surprising direction for G to act in is due to the definition Eq. (2.5) below (cf. 

the relation between the active and passive view of rotations). In order to give a general 
form of the effective low-energy Hamiltonian, we define another operator X{n, /j,, A), 

X{n, n, A) = (1 + G) (1 + G^G + GG^)^. (2.3) 

{n is a real number.) The inverse of X{n, /i, A) can be obtained explicitly, 

X-\n, n, A) = (1 + G^G + GG^y{l - G). (2.4) 
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The special case n = is equivalent to the transformation introduced in Ref. [24] to 
relate the hermitian and non-hermitian effective operators in the energy-independent 
Suzuki-Lee approach. We now consider the transformation of H{A) defined as 

H{n, n, A) = X-\n, A)H{A)X{n, //, A), (2.5) 

where we have 

H{A) H{n, 12, A) = H{i2) + 6H{i2, A). (2.6) 
One can prove that if H{n, /i, A) satisfies the desirable decoupling property, 

QH{n,ii,A)P^O, (2.7) 

or more explicitly, by substituting the definition of X{n,ii,A) and X~^{n, /i, A) from 
Eqs. (2.3-2.4), 

QH{A)P + QH{A)QG - GPH{A)P - GPH{A)QG = 0, (2.8) 

that H^^[ii) = TC{n, ji) = PH{n, ji, A)P is an effective Hamiltonian for the low energy 
degrees of freedom. In other words, it should have the same low-energy eigenvalues as 
the original Hamiltonian. The proof is as follows: 

Consider an eigenvalue equation in the P space with {\(j){k)) e P}, 

PH{n, II, A)P\(l)(k)) = EkPX-\n, ii, A)X{n, A)P\(t){k)). (2.9) 

By multiplying both sides by X{n,ii,A) and making use of the decoupling property 
Eq. (2.7), we obtain 

H{A)X{n, li, A)P\cf)(k)) = EkX{n, A)P|0(A;)). (2.10) 

This equation means that in Eq. (2.9) agrees with one of the eigenvalue of H[A) 
and X{n, /i, A)P\4>{k)) is the corresponding eigenstate. Now we demand that 

H'''{fi) = H{fi) + Hc{fi). (2.11) 

This requirement uniquely determines the counterterm He- 

In the same way we can also obtain the Q-space effective Hamiltonian, from the 
definition of H{n, fi. A). It can be seen that if G satisfies the requirement in Eq. (2.8), 
then we have additional decoupling condition 

PH{n,ii,A)Q ^0. (2.12) 
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We will argue later that both of the decoupling conditions Eqs. (2.7) and (2.12) are 
necessary in order to have a sector-independent renormalization scheme. The word 
"sector" here means the given truncated Fock space. Let us now clarify the meaning 
of this concept. To maintain the generality of the previous discussion, we use here the 
well known Bloch-Feshbach formalism [12, 25, 26]. The Bloch-Feshbach method exploits 
projection operators in the Hilbert space in order to determine effective operators in 
some restricted model space. This technique seems to be more universal than Wilson's 
renormalization formulated in a Lagrangian framework. This is due to the fact that in 
the Bloch-Feshbach formalism, other irrelevant degrees of freedom (such as high angular 
momentum, spin degrees of freedom, number of particles, etc.) can be systematically 
eliminated in the same fashion. 

Assume that the full space Schrodinger equation is if(A)|'0) = Elip) and for sim- 
plicity l'^) has been normalized to one. We explicitly construct the energy dependent 
effective Hamiltonian in this formalism, 

^eff ^ p-Jjp ^ p-JjQ \- QHP, (2.13) 

where H can be a similarity transformed Hamiltonian. This equation resembles the 
one for Brueckncr's reaction matrix (or "G" -matrix) equation in nuclear many-body 
theory (NMT). In the same way for arbitrary operator O (after a potential similarity 
transformation), we construct the effective operator 

0- - POP + PHQ^J^QOP + POQ^-L^QHP 

The ^'-dependence in Eqs. (2.13) and (2.14) emerges from the fact that the effective 
interaction in the reduced space is not assumed to be decoupled from the excluded space. 
However, by using the decoupling conditions introduced in Eqs. (2.7) and (2.12), we 
observe that energy dependence can be removed, and the effective operators become 

H^^ = PHP ^n{n,fi), 

O^^ ^ POP^O{n,n). (2.15) 

The decoupling property makes the operators in one sector independent of the other sec- 
tor. The effects of the excluded sector is taken into account by imposing the decouphng 
conditions. This is closely related to the folded diagram method in NMT for removing 
energy-dependence [27]. (It is well-known in NMT that ^'-dependence in the G-matrix 
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emerges from non-folded diagrams which can be systematically eliminated using the 
effective interaction approach). The above argument was given without assuming an 
explicit form for X and thus the decoupling conditions are more fundamental than the 
prescription used to derive these conditions. We will show later that one can choose 
a transformation X, together with the model space and its complement, which avoids 
the occurrence of "small energy denominators" . We now show that Lorentz covariance 
in a given sector does not hinge on a special form of similarity transformation. We 
assume ten Poincare generators Lj satisfying 

[L,,L^]^Y.'''iLk, (2.16) 

where the a^j are the known structure coefficients. One can show that if the operators 
Li satisfy the decoupling conditions QL^P — and PL^Q — then it follows that 

[Lf,Lf]=J24Lf. (2.17) 

This leads to a relativistic description even after simultaneously integrating out the 
high-frequency modes and reducing the number of particles. Indeed we conjecture that 
requiring the decoupling conditions makes the effective Hamiltonian free of Lorentz- 
noninvariant operators for a given truncated sector regardless of the regularization 
scheme. However, one may still be faced with an effective Hamiltonian which violates 
gauge invariance (for e.g., when sharp cutoff is employed). 

Note that the solution to Eq. (2.8) is independent of the number n. One can make 
use of Eq. (2.8) and its complex conjugate to show that for any real number n, the 
following relation for the effective low-energy Hamiltonian 

n{n,n) =n\-n-l,fi). (2.18) 

The case n = — 1/2 is special since the effective Hamiltonian is hcrmitian, 

n(-l/2,ii) = (P + G^Gy/^H(A)(P + G)(P + G^G)-^/^. (2.19) 

Hermiticity can be verified from the relation [28] 

e^P = {l + G){P + G^G)-^/\ (2.20) 

where, 

T = arctan(G' - G^) = J] ^--^(G(G'^G')" - h.c). (2.21) 

n=0 

Since the operator T is anti-hermitian, is a unitary operator. Prom the above 
expression Eq. (2.19) can be written in the explicitly hermitian form 

n(-l/2,ii)^ Pe-^H{A)e^P. (2.22) 
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As was already emphasized, renormalization based on unitary transformations is more 
complicated and non-economical. Thus we will explore a non-unitary approach. An 
interesting non-hermitian effective low-energy Hamiltonian can be obtained for n = 0, 

niO,ii,A)^PH{A){P + QG). (2.23) 

This form resembles the Bloch and Horowitz type of effective Hamiltonian as used in 
NMT [25], and was used by Wilson in quantum field theory [12]. In the context of 
the CCM, this form leads to the folded diagram expansion well known in many-body 
theory [29]. It is of interest that various effective low-energy Hamiltonians can be 
constructed according to Eq. (2.15) by the use of the mapping operator G which all 
obey the decoupling property Eq. (2.7) and Eq. (2.12). Neither perturbation theory 
nor hermiticity is essential for this large class of effective Hamiltonians. 

The CCM approach is, of course, just one of the ways of describing the relevant 
spectrum by means of non-unitary transformations. According to our prescription the 
model space is P : {\L) ^\0,b)h, L < /i}, where \0,b)h is a bare high energy vacuum 
(the ground state of high-momentum of frce-Hamiltonian) which is annihilated by all 
the high frequency annihilation operators {C/} (for a given quantization scheme, e.g., 
equal time or light-cone) , the set of indices {/} therefore defines a subsystem, or 
cluster, within the full system of a given configuration. The actual choice depends 
upon the particular system under consideration. In general, the operators {C/} will 
be products (or sums of products) of suitable single-particle operators. We assume 
that the annihilation and its corresponding creation {C|} subalgebras and the state 
\0,b)h are cyclic, so that the linear combination of state {C]|0,&)^} and { h{b,0\Ci} 
span the decoupled Hilbcrt space of the high-momentum modes, {\H)}, where fi < 
H < A. It is also convenient, but not necessary, to impose the orthogonality condition, 
(0|C/Cj|0) = 6{I,J), where 6{I,J) is a Kroncckcr delta symbol. The complement 
space is Q : {\L) ^{\H) — \0,b)h)}- Our main goal is to decouple the P-spacc from 
the Q-space. This gives sense to the partial diagonaUzation of the high-energy part of 
the Hamiltonian. The states in full Hilbert space are constructed by adding correlated 
clusters of high-energy modes onto the P-space, or equivalently integrating out the 
high-energy modes from the Hamiltonian, 

I/) =X(/.,A)|p) =eV^'|0,6),(g)|L) =e^|0,6),(g)|L), (2.24) 
{f\ = (^1 (8) h{b, 0\X-\^i, A) = (L| (g) H{0\e''e-', (2.25) 

where the operators X{ii,A) and X~^{ii,A) have been expanded in terms of indepen- 
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dent coupled cluster excitations / 



m 




m=0 



) 



m 



(2.26) 



m=0 



I 



Here the primed sum means that at least one fast particle is created or destroyed 



generally commutable in the low-energy Fock space, whereas they are by construction 
commutable in the high-energy Fock space. It is immediately clear that states in the 
interacting Hilbert space are normalized, (/|/) = 0|0, = 1. We have two types 
of parameters in this procedure: One is the coupling constant of the theory (A) , and 
the other is the ratio of cutoffs (/x/A). The explicit power counting makes the degree 
of divergence of each order smaller than the previous one. According to our logic, 
Eq. (2.23) can be written as 



with X(/x,A) and X'^{^ji,K) defined in Eq. (2.24) and Eq. (2.25). We require that 
effective Hamiltonian 7i(/x) obtained in this way remains form invariant or coherent [10, 
23]. This requirement satisfies on an infinitely long renormalization group trajectory 
and thus does constitutes a renormalized Hamiltonian. Thereby one can readily identify 
the counter terms produced from expansion of Eq. (2.27). 

It is well-know in many-body applications that the exponential Ansatz Eqs. (2.24) 
and (2.25) guarantees automatically proper size-exclusivity and conformity with the 
Goldstonc linked- cluster theorem to all level of truncation. This paramctrization 
does not manifestly preserve hcrmitian conjugacy. However, it is compatible with the 
Hellmann-Feynman theorem (HFT), in other words, demanding hermiticity will violate 
this theorem at any level of truncation. On the other hand, with this parametrization 
the phase space {s/, s'j} for a given m is a symplectic differentiable manifold. Thereby 
all the geometrical properties of the configuration space can be precisely defined. (In a 
more ordinary language, the canonical equations of motion with respect to phase space 
define a set of trajectories, which fill the whole dynamically allowed region of the phase 
space.) There is a deep connection between these three properties and we can not give 
up one without losing something else as well [30]. The individual amplitudes for a given 
m, {sf, s'™} = {si, S/}m, can generally be functionals of the low- and high-energy field 
operators and have to be fixed by the dynamics of quantum system. This is a com- 
plicated set of requirements. However, we require less than that. Suppose that after 



(/ 7^ 0), and momentum conservation in P Q is included in si and s'j. S, 



m 



{S'J are not 



nifx) = h{b,0\X-\f,,A)H{A)X{f,,A)\Q,b)h, 



(2.27) 
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a similar transformation of Hamiltonian, H, we obtain an effective Hamiltonian of the 
form 

H = //(low) + i/free(high) + CjVijCj, (2.28) 

where Vjj is an arbitrary operator in the low frequency space. The / and J in- 
dices should be chosen such that the last term in Eq. (2.28) contains at least one 
creation- operator and one annihilation-operator of high frequency. By using Rayleigh- 
Schrodinger perturbation theory, it can be shown that the free high-energy vacuum 
state of i7free(high) is annihilated by Eq. (2.28) and remains without correction at any 
order of perturbation theory. Having said that, we will now consider how to find the 
individual amplitudes {s/, s'j}„i that transfer the Hamiltonian into the form Eq. (2.28). 
We split the Hamiltonian in five parts: 

H^H, + //^■^^^(high) + VciCj) + Va{Ci) + Vb, (2.29) 

where Hi contains only the low frequency modes with k < fi, H2 is the free Hamiltonian 
for all modes with fi < k < A, Vc contains low frequency operators and products of 
the high frequency creation operators C] and Va is the hermitian conjugate of Vc- The 
remaining terms are contained in Vb, these terms contain at least one annihilation and 
creation operators of the high energy modes. Our goal is to eliminate Vc and Va since 
Vb annihilates the vacuum. The ket-state coefficients {si}m are worked out via the 
ket-state Schrodinger equation H{A)\f) — E\f) written in the form 

(0|C^e"%e^|0) = 0, V/ 7^ 0. (2.30) 

The bra-state coefficients {§1, s'j}jn arc obtained by making use of the Schrodinger 
equation defined for the bra-state, {f\H{A) = {f\E. First we project both sides on 
C||0), then we eliminate E by making use of the ket-state equation projection with the 
state {0\e^'Cl to yield the equations 

{0\e^'e-^[H, Cj]e^e-^'|0) = 0, V7 ^ 0. (2.31) 

Alternatively one can in a unified way apply e^e~^' Cl\0) on the Schrodinger equation 
for the bra-state and obtain 

(0|e^'e-%eV^'cJ|0) = 0, V7 ^ 0. (2.32) 

Equation (2.30) and Eqs. (2.31) or (2.32) provide two sets of formally exact, micro- 
scopic, operatorial coupled non-linear equations for the ket and bra. One can solve the 
coupled equations in Eq. (2.30) to work out {si}m and then use them as an input in 
Eqs. (2.31) or (2.32). 
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It is important to notice that Eqs. (2.30) and (2.31) can be also derived by re- 
quiring that the effective low-energy Hamiltonian defined in Eq. (2.27), be station- 
ary (i.e. 51-L{ii) = 0) with respect to all variations in each of the independent func- 
tional {si,s'j}rn- One can easily verify that the requirements 5Tl{iJ,)/5si — and 
5n{n)/5s'j = yield Eqs. (2.30) and (2.31). The combination of Eqs. (2.30) and (2.31) 
does not manifestly satisfy the decoupling property as set out in Eqs. (2.7) and (2.12). 
On the other hand Eqs. (2.30) and (2.32) satisfy these conditions. Equations (2.30) and 
(2.32) imply that all interactions including creation and annihilation of fast particles 
("/") are eliminated from the transformed Hamiltonian 7i(/i) in Eq. (2.27). In other 
words, these are decoupling conditions leading to the elimination of Vc and Va from 
Eq. (2.29), which is, in essence, a block-diagonahzation. Therefore it makes sense for 
our purpose to use Eqs. (2.30) and (2.32) for obtaining the unknown coefficients, losing 
some of the elegance of the CCM elsewhere. 

So far everything has been introduced rigorously without invoking any approxima- 
tion. In practice one needs to truncate both sets of coefficients {sj, s'j}m at a given 
order of m. A consistent truncation scheme is the so-called SUB(n, m) scheme, where 
the n-body partition of the operator {S, S'} is truncated so that one sets the higher 
partition with / > n to zero at a given accuracy m. Notice that, Eqs. (2.31) and (2.32) 
provide two equivalent sets of equations in the exact form, however after the trunca- 
tion they can in principle be different. Eqs. (2.30) and (2.32) are compatible with the 
decoupling property at any level of the truncation, whereas Eqs. (2.30) combined with 
(2.31) are fully consistent with HFT at any level of truncation. Thus the low-energy 
effective form of an arbitrary operator can be computed according to Eq. (2.15) in the 
same truncation scheme used for the rcnormalization of the Hamiltonian. In particular, 
we will show that only in the lowest order (m = 0), equations (2.31) and (2.32) are 
equivalent, independent of the physical system and the truncation scheme. 

Although our method is non-perturbative, perturbation theory can be recovered 
from it. In this way, its simple structure for loop expansion will be obvious and we 
will observe that at lower order hermiticity is preserved. Now we illustrate how this is 
realizable in our approach. Assume that Vc and Va are of order A, we will diagonalize 
the Hamiltonian, at leading order in A up to the desired accuracy in /x/A. We use the 
commutator-expansion 

e-'He' ^H + [H,S] + ^[[H,S],S] + ... . (2.33) 

Eq. (2.30) can be organized perturbatively in order of m, aiming at elimination of the 
high momenta degree of freedom up to the first order in the coupling constant, thus 



- 11 - 



yields 

m^0:{0\Ci(Vc + [H2,So])\0)^0, 

m = 1 : {0\Cj([H,, So] + [H2, Si] + [Va, Si] + [Vc, Si])\0) = 0, 

m = n : (0|C,([ifi, ^-^-i] + [H2, Sn] + [Va, Sn] + [Vc, |0) = 0, (2.34) 

where J 7^ 0. Notice that 5*0 is chosen to cancel Vc in the effective Hamiltonian, hence 
it is at least of order of A, consequently it generates a new term [Hi, So] which is of 
higher order in ^/A and can be canceled out on the next orders by Si. The logic for 
obtaining the equations above is based on the fact that Sn should be smaller than 
^^-i (for sake of convergence) and that the equations should be consistent with each 
other. Since H2, Va, Vc ~ A and Hi pa /x, from Eq. (2.34) we have the desired relation 
Sn ~ The same procedure can be applied for Eq. (2.32) which leads to the 

introduction of a new series of equations in order of m, 

m = 0:{0\{VA-[H2,S'o])C]\0) = 0, 

m - 1 : mHi, S',] + [H2, S[] + [Vc, S[] + [Va, S[] - [Va, Si])C]\0) = 0, 

m^n: (0|([i/i, S'^_,] + [H^, S'^] + [Vc, S'J + [Va, S'^] - [Va, S^])Cj\0) = 0(2.35) 

Alternatively, we can use Eq. (2.31) to yield the equations 

m = : (0| {[Va, CJ] - [[H^, Cj], S',] ) |0) = 0, 

m = l:(0|([[y^,Cj],^i] - [[H2,Cl],S[] - [[Va,C\], S[])\0) ^ 0, 

m = n:{0\{[VA,Cj],Sn] - [[H2,Cj], S'^] - [[^4, C]], 5^] ) |0) = 0. (2.36) 

It is obvious that at order m — 0, Eqs. (2.35) and (2.36) are the same and S'q — S^, 
which indicates that the similarity transformation at this level remains unitary. It 
should be noted that diagonalization at first order in the coupling constant introduces 
a low-energy effective Hamiltonian in Eq. (2.27) which is valid up to the order A^. In 
the same way, diagonalization at second order in A modifies the Hamiltonian at order 
A"^ and leads generally to a non-unitarity transformation. In this way one can proceed 
to diagonalize the Hamiltonian at a given order in A with desired accuracy in /i/A 
. Finally, the renormalization process is completed by introducing the correct Z{A) 
factors which redefine the divergences emerging from Eq. (2.27). 
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3. Conclusion and outlook 

In this paper we have outhned a strategy to derive effective renormahzed operators in 
the Hamiltonian formulation of a field theory. An example of the application of the 
method can be found in Ref. [21]. The effective low-frequency operator is obtained by 
the condition that it should exhibit decoupling between the low- and high-frequency 
degrees of freedom. We showed that the similarity transformation approach to renor- 
malization can be systematically classified. The non-hermitian formulation gives a very 
simple description of decoupling, leading to a partial diagonalization of the high-energy 
part. In the recent paper [21], we showed that non-unitarity representation inherent 
in our formulation is in favour of economic computation and does not produce any 
non-hermiticity in the relevant terms. The techniques proposed are known from the 
coupled cluster many-body theory and invoke neither perturbation nor unitarity trans- 
formation. We showed that our formalism can be solved perturbatively. In this way, it 
was revealed that diagonalization at first order in coupling constant defines a correct 
low-energy effective Hamiltonian which is valid up to the order A^. One can show that 
the non-hermiticity of the effective Hamiltonian is controllable and might appear in 
higher order which is beyond our approximation or in irrelevant terms which can be 
ignored in renormalization group sense. 

One of the key features which has not yet been exploited is the non-perturbative 
aspect of the method; it may well be able to obtain effective degrees of freedom that 
are very different from the ones occurs at the high-energy scale. This is a promising 
avenue for future work. 
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